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Abstract 

We show that the Kazakov-Migdal (K-M) induced gauge model in d di- 
mensions describes the high temperature hmit of ordinary lattice gauge theo- 
ries in d-l- 1 dimensions. The matter fields are related to the Polyakov loops, 
while the spatial gauge variables become the gauge fields of the K-M model. 
This interpretation of the K-M model is in agreement with some recent results 
in high temperature lattice QCD. 
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1. Introduction 



Recently V.Kazakov and A.Migdal proposed a new lattice gauge model, in which 
the gauge self-interaction is induced by scalar fields in the adjoint representation 
The action they propose is defined on a generic d- dimensional lattice and has the 
following form: 

S = J2NTT[m^(j)'^{x) -J2(Pix)U{x,x + ^i)(j){x + ^m)U^{x,x + iJ,)] (1) 

where is an Hermitian N x N matrix defined on the sites x of the lattice and 
U{x,x + /i) is a unitary N x N matrix, defined on the links {x,x + /i), and plays 
the role, as in the usual lattice discretization of Yang-Mills theories, of the gauge 
field. Integrating over the scalar field one can induce an effective action for the 
gauge field, 

J DUD^exp{-S) ~ J DUexp{-Sind[U]) (2) 

with: 

s \u]- (^) 

^rad[U\ 2 V^[r](2m2)nr] ' ^"^^ 

where /[F] is the length of the loop F, [/[F] is the ordered product of link matrices 
along F and the summation is over all closed loops. 

The model turns out to be solvable in the large limit, for any space-time 
dimension d [0, 0, |^ and exactly solvable, for any value of in d=l IQ. A further 
interesting feature of the model is its deep connection with the theory of non- 
critical strings, in particular it can be shown that in the d = 1 case the K-M 
model is equivalent to the vortex free sector of the d=l string Despite all these 
nice features the model seems to miss the original goal of Kazakov and Migdal, 
since the induced gauge theory (^, due to the fact that the matter fields are in 
the adjoint representation, has a super-confining behaviour and seems to have 
the wrong perturbative vacuum 0, ^. Several improvement of the original K-M 
action have been proposed to avoid these problems ||, |, 0, |Tl|] , but even if many 
interesting results have been obtained the problem of the identification with a pure 
lattice gauge theory of the ordinary type (namely with the ordinary confinement 
behaviour) while keeping exact solvability is still open. 

In this letter we suggest a new point of view on the problem and show that 
the K-M model can be identified before integrating over the matter fields with an 
ordinary SU{N) lattice gauge theory in the high temperature limit, the matter 
fields being related to the Polyakov loops, while the spatial gauge variables become 
the gauge fields of the K-M model. 

This letter is organized as follows: after a short introduction (just to set nota- 
tions) on finite temperature lattice gauge theory (sect. 2), we describe and discuss 
the equivalence with the K-M action in sect. 3 and 4. Sect. 5 is devoted to the 
continuum limit and sect. 6 to some concluding remarks. 
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2. Finite temperature lattice gauge theories 

Let us consider a pure gauge theory with gauge group SU (N) defined on a + 1 
dimensional cubic lattice. In order to describe a finite temperature LGT we must 
take periodic boundary conditions in one direction (which we shall call from now 
one "time-like" direction), while the boundary conditions in the other d direction 
(which we shall call "space-like") can be chosen freely. Let us take a lattice of Nt 
(Ng) spacings in the time (space) direction. The theory will contain only gauge 
fields described by the link variables Un;i G SU{N) where n = {x, t) denotes the 
space-time position of the link and i its direction. It is useful to choose different 
couplings in the time and space directions. Let us call them and (3s respectively. 
Let us take the simplest choice for the lattice gauge action, namely the Wilson 
action: 

Sw = Y. N RelptY. tr{Un■,o^) + /5s E tr{Un;,) \ (4) 

n [ i i<j J 

where Un-oi {Un;ij) are the time-like (space-like) plaquette variables, defined as 
usual: 

Un,, = U^,JJ^^,pl^.Ui, (5) 

Ps and Pt are related to the (bare) coupling constant g and the temperature T 
of the gauge theory by the relations 0: 

2 _.3-d 



where a is the space-like lattice spacing while is the time-like spacing. The two 
are related by the adimensional ratio e = j^^- We can solve the above equations 
in terms of e as follows: 

A = ^a'-' (7) 

re 



2e 
? 

In a finite temperature discretization it is possible to define gauge invariant vari- 
ables which are topologically non-trivial loops, closed due to the periodic boundary 
conditions in the time directions. The simplest choice is the Polyakov loop defined 
as follows: 

Nt 

P{x)=trX{{U^^t,o) (9) 

t=i 



^ Notice that the couphng constant g is rescaled here by a factor N with respect for instance 
to the notations of ref. [fl3| 
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where x labels the space coordinates of the lattice sites. Moreover, an important 
feature of the finite temperature theory with respect to the zero temperature case 
is that it has a new global symmetry (independent from the gauge symmetry) with 
symmetry group the center C of the gauge group (in our case Z]^). The Polyakov 
loop turns out to be a natural order parameter for this symmetry. 

It is possible to obtain, just from the definition itself of the model, some general 
properties in the high temperature regime (see for instance In this region the 

symmetry with respect of the center of the group is broken, the theory is deconfined, 
the Polyakov loop has a non-zero expectation value and, what is more important, 
it is an element of the center of the gauge group. Physically this means that the 
links in the time direction fluctuate around one of the N^^ roots of unity of the 
Zjq group. We can lift this degeneracy by adding to the Wilson action a "magnetic 
term" 

Nt 

Sm = Nh J2Retrl[{U,,t,o) (10) 

X t=l 

and eventually send h 0. This procedure selects the vacuum in which 
< P{x) >= 1. 

The spatial links are static up to gauge transformations and this fact tells us 
that the spatial degrees of freedom of the model behave as if they belonged to a zero 
temperature model in one dimension less and with coupling constant g^T. If one 
could integrate out these spatial degrees of freedom the resulting effective action for 
the Polyakov loops would be that of a d- dimensional spin model with short range 



spin-spin interaction 12 



3. The K-M model as an high temperature LGT. 

Let us take the d + 1 dimensional LGT described above (with gauge group 
SU{N)) and consider first a very special case, namely A^t = 1. We will be interested 
in the T — oo (hence e — 0) limit. 

Let us call St {Ss) the contribution of the temporal (spatial) plaquettes to the 
action, with a "magnetic term" of the type (p!0| ) included in St. Due to the boundary 
conditions and to the Nt = I position we can rewrite St as: 

St = N(3tY: tr {Ux,^U,^i,oUlMo} + NhY.Re {U,,o) (11) 

x,i X 

This model is interesting in itself and represents the generalization of the K-M 
model to "matter" fields described by unitary matrices. Notice also that since the 
space contribution Ss to the action is depressed as ~ 1/T, in the T ^ oo limit 
this model should capture most of the features of high temperature LGT. 
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Following Svetitsky and Yaffe one can argue at this point that in the high 
temperature limit each temporal plaquette fluctuates around the identity . It is 

easy to see that the fluctuations are of order • In fact if we put Upi = e ^ 

with (j) a traceless hermitian matrix we have: 

dUpie""^''' = / c/0e-^ '''' [I + 0(1/A)] (12) 



Correspondingly the link variables in the time direction fluctuate around one of 
the A^*'' roots of unity of the Z^r group with fluctuations still of order I/a/S*. When 
the magnetic term is switched on the vacuum is forced to be U^-o = 1 and we can 
expand U^-q in the following way: 

i^a (j){x) (t?{x) , , 

CI,,.ev^=l + .^-^ + ... (13) 
By inserting ([T3|) in (|l^) we find: 



St = f^ ^ tr id Q^l + j^Y. (-m^<P{xf + ^ U^.,(P{x)Uy{x + l)\ | (14) 



where ~ ^ + 2I1 volume of the d dimensional space. In ref. [|r2| only 

the constant leading contribution of order /3j was considered; here we recognize that 
the contributions of the fluctuations, of order 1, coincide, apart from an irrelevant 
constant, with the K-M action with arbitrary (eq. (P) and in the limit of zero 
magnetic field with the K-M action at the critical point = d. The expansion 
(p^) and in particular the choice of the expansion parameter — ^ is also justified 



a posteriori by the analysis of the distribution of the eigenvalues of 0(x) resulting 
from the action {\L^ . It was shown in |^ that at least in the large limit and for 
> d the eigenvalue distribution of (j){x) is semicircular with a finite , non-zero 
radius. This implies that the eigenvalues of the original unitary matrix U^-q are 
restricted to a region of order around 1 on the unit circle, thus justifying the 



expansion (p!3|). For d > I such distribution of the eigenvalues survives also in the 
limit h (with h > 0) namely d, whereas for d = 1 the radius of the 

eigenvalue distribution goes to infinity as h —>■ (m^ ^ 1) • This shows that in the 
one dimensional case the critical point can only be described in the unitary matrix 
formulation given by (0). 

Let us consider now the case of an arbitrary Nt. The action (without the con- 
tribution of spatial plaquettes) is: 

Nt 

St = iVA E E tr{U,,;JJ,^,;oUh+iMAo} + NhJ2 ReP{x) (15) 

x,i t=l X 
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where P{x) is the Polyakov loop given in (|^). It is possible to choose the gauge 
in such a way that Ug^t;0 = 1 for ^ = 1)2, ..,Nt — 1 so that the only non trivial 
links in the time direction are Ug^Nt;0 ? related to the Polyakov loop by the relation 
P{x) = trUni^Nffi- The plaquettes at t < Nf then reduce to U^^t;iU'^gt+i-i ■ follows 
that in the large (3t limit the spatial link variables at different values of t coincides 
up to fluctuations of order So we can put 



(16) 



where we are free to choose tpi 
we also have: 

U. 



X, 1) = 0. In complete analogy with the case A''^ = 1 



exp{i 




(17) 



where the factor y/Nt in the exponent at the r.h.s. of ( p!7|) is needed , as we shall 
see, for a smooth Nt ^ oo limit. By inserting equations (0) and (p!7| ) into the 



action ([T5|) and by expanding in powers of l//3f up to terms of order 1 , one obtains 
an action quadratic in the ipii^yVY^'- 



St = Ntr J2{Nt(3t + E IM^^ t+1)- ^Ux, t)] - ^K^, Nt) 



t=2 



Ntij,{x,Nt)[(f){x + i) 
Nt 



m 



(X + i] 



-uyix)u,. 



'-Nh-^ct>\x)} 

Pt 



:i8) 



The fluctuations tpi{x,t) can be eliminated by performing the corresponding gaus- 
sian integrals and the flnal result coincides with the K-M model as given by eq.(|T^[|. 
The same result can be obtained also by noticing that, as the self interaction term 
is neglected, each spatial link variable U:g^t;i belongs to just two temporal plaquettes 
and that the corresponding integral can be done by using character expansion in 
analogy to lattice QCD in two dimensions. This allows to deflne a "renormalized" 
coupling Peff through the relation: 



eff) 



Nt 



(19) 



where If and Jq are the coefficients of the fundamental and identity characters in the 
character expansion of the action. Then it can be shown that the large (3 behaviour 
of (5eff{P) is as expected [3eff 



Nt 



m) 



(see for instance ref. 
The final result eq.(0) is then independent from the number Nt of lattice spacing 
in the time direction provided the field (f){x) is normalized as in (|l^. Notice that 



^Notice that fix in ([T^ ) is now replaced by the volume of the whole space-time lattice 
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the r.h.s. of ( |T^) is also independent of Nt according to the definition of Pt given 
in (^. So the continuum hmit in the time dimension A^j — oo is well defined and 
leads to the K-M model. 

Hence we can conclude that 

for d > 1 the d-dimensional K-M model is a good description of the small 
fluctuations of the Polyakov loops around their minimum (frozen) position in the 
high temperature limit of a pure lattice gauge theory in d + 1 dimensions with the 
Wilson action. 



4. Comments and remarks 

The integration over the fields in (|1^) leads to the induced gauge action eq. 
which can now be interpreted as describing the spatial degrees of freedom of a pure 
LGT at T = oo. Since, as mentioned above, the space degrees of freedom in the 
high temperature limit behave as a zero temperature d-dimensional gauge theory 
with coupling constant g'^j-jr = g^T, the K-M model corresponds to the geff = oo 
point of the theory thus leading to the superconfinement behaviour. 

At T < oo one has to add a small (order 1/T) gauge self-interaction term. This 
term obviously destroys the exact solvability of the model, but it can be treated as 
a small perturbation around the K-M solution. The main outcome of our analysis 
is probably in the fact that we now understand that such a perturbative analysis 
would be an high-temperature expansion for the LGT. 

Let us consider for instance the expectation value of a non backtracking Wilson 
loop which is zero at T = cxd due to the superconfinement. Its first non trivial 
perturbative contribution in the 1/T expansion is obtained by filling the loop by 
elementary plaquettes 0. The resulting "filled Wilson loops" is invariant under the 
local Zn symmetry and has to be evaluated using the measure of the pure M-K 
induced action. This was done in ref. [0 where it was shown that the "filled Wilson 
loop" has an area law behaviour. It is also remarkable and not quite understood 
yet that in this context intriguing connections with 2d spin models emerge. Our 
interpretation is that the "filled Wilson loop " is the first non vanishing term in 
a large temperature expansion of the vacuum expectation value of the ordinary 
Wilson loop. 

The K-M model is only apt to describe fluctuations around one given vacuum of 
the high temperature LGT; in other words it can only describe finite temperature 
LGT in the broken phase. However, even in the framework of the K-M model 

^In more than two spatial dimension one has to sum over the contributions of aU surfaces made 
of elementary plaquettes and having the original loop as a boundary. These contributions are of 
order T~-^ where A is the area of the surface, so that the surface of minimal area dominates at 
high temperature. 
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one has at least a qualitative understanding of the restoring of the Zjv symmetry 
below the critical temperature. In fact the expectation value of the Polyakov loop 



< p(f) >=< Tre V > (20) 

can be evaluated in the large N limit by replacing 0(x) with its classical value given 
in ref. 0. At the critical point = d this is given by a semicircular distribution 
of eigenvalues of radius r = y^^^Eij- At high temperature, more precisely for 

^ TT, the eigenvalues are peaked around one vacuum and the symmetry 

is broken. However for values of (3t such that rJ ^ ^ tc the eigenvalues of J ^(j){x) 



are distributed over several vacua on most of the unit circle and so the symmetry 
is eventually restored. 

Without a more detailed understanding of the symmetry restoration nothing 
new can be said in this context on the Svetitsky-Yaffe conjecture that the critical 
behavior of a ((i+l)-dimensional finite temperature gauge theory is the same as that 
of a (i-dimensional spin model with the same Z^r symmetry. At this stage at least 
the emergence of spin models 0] in the expectation value of the Wilson loop does 
not appear to be correlated with the Svetitski-Yaffe conjecture. It should be noticed 
however that if the contribution of the spatial plaquettes is neglected, an effective 
theory for the Polyakov loop valid also in the region of large (f){x) could in principle 
be obtained directly from ( |T3| ) by integrating over the spatial links variables. This 
would require the use of a generalization of the Itzykson-Zuber formula described 
in 1^ which is exact for N = 2 and valid as an asymptotic formula for N > 2. 

The instability which occurs in the < d region, corresponding to negative 
values of h in ([ISD , can be understood from the high temperature LGT point of 
view as a consequence of the fact that with our "magnetic field" in this regime we 
are actually pushing the system out of the chosen vacuum which becomes unstable. 
Similarly the instability due to the presence of a linear term ||I4| which occurs if 
one looks at U(N) instead of SU(N) K-M models, can be understood as a signature 
of the Goldstone modes of the broken U(l) symmetry which one has in this case. 

Let us finally consider the case d = 1. Spatial plaquettes are absent in this 
case and eq. (JT^ gives the complete action irrespective of the value of jSt- This is 
just two dimensional QCD on a cylinder or , in case the spatial dimension is also 
compactified, on a torus. We know from ref. that in the d = 1 case at the critical 
point rn? = 1 the distribution of the eigenvalues of Uft^Nffl is not confined to a small 
region of the unit circle. As a consequence we are always in the unbroken phase 
and the description of the system in terms of the K-M model is not valid. On the 
other hand by expanding both spatial and temporal links in a fashion similar to 
eq.(]T6D one can reduce the action of eq. ( |T5| ) with periodic boundary conditions in 
both space and time to an action containing only one space and one time-like link. 



in agreement with general result on lattice QCD in two dimensions [ITT 
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5. Taking the continuum limit 



In the previous sections we have derived the K-M model as a high temperature 
limit of lattice QCD. However while in ordinary lattice QCD one has a well defined 
continuum limit, this is not the situation in the K-M model with a quadratic poten- 
tial at the critical point = d. It was pointed out in p that if one approaches 

= d from below there exists a solution of the master field equation correspond- 
ing to a semicircular distribution of eigenvalues whose radius r goes to infinity as 

— > (i, a signal of the existence of a critical point. Such configuration however is 
a local maximum of the free energy; moreover the free energy itself is unbounded 
from below for < d. This can be cured by adding a higher order term to the 
action, for instance a quartic term jTrcj)^. From the point of view of the high 
temperature expansion this is not such an ad hoc adjustment as it might appear at 
first, in fact the continuum limit corresponds to the limit where the radius of the 
eigenvalue distribution goes to infinity and in that limit higher order terms in 0(x) 
are expected to be relevant. Therefore in the continuum limit the high temperature 
QCD is most likely described by a K-M model with a higher order potential in 



agreement, as discussed below, with some recent results on lattice QCD [0, |T7 
It has indeed been shown in that, with the addition of a quartic term ^Trcp^ 
to the action, a second order critical point can be reached, at least in the case of 
SU(2) in d=4. The analysis of the phase diagram can be done rather easily within 
a mean field approximation: in the A, plane one finds a line of first order tran- 
sitions starting from the K-M critical point and ending with a second order critical 
point located at A = 2.57 and = 2.260. Moreover direct Montecarlo simulations 
show that these mean field results are quite accurate and even the critical point 
location is essentially confirmed by the computer simulation. We have made a sim- 
ilar mean field analysis in the d = 3 case, in which we are interested, showing a 
similar scenario: a line of first order phase transitions starting from = 3 and 
ending with a second order critical point located at A ~ 2.2, ~ 1.56. In ^5 



it was stressed that the corresponding continuum theory had nothing to do with 
ordinary SU(2) gauge theory. Our suggestion is that it should instead describe the 
high temperature, deconhned phase of SU(2) (with one more space-time dimen- 
sion). This conjecture is in remarkable agreement with some recent results on the 
high temperature behaviour of lattice gauge theories, both within a perturbative 
approach |T^, and with montecarlo simulations [|1^. The main idea behind |16, |1^ 



is that (3+1) dimensional gauge theories at high temperature undergo a peculiar 
form of dimensional reduction. The original picture [|18[ of a complete dimensional 
reduction (namely a complete decoupling of the degrees of freedom in the compacti- 
fied time direction, which would lead to an effective three-dimensional theory of the 
ordinary type) does not occur in general and one finds instead a three dimensional 
gauge theory coupled with a scalar field (p in the adjoint representation with a non 

^Notice that there is a factor of two between our normaUzation of the mass parameter and that 

of 
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trivial potential V{(f)). This potential has a quadratic term with the wrong sign to 
be identified as a mass term and a quartic contribution which stabilizes the action, 
exactly as in our mean field solution. Moreover the action used in (eq.(20) 
of fl^) for the montecarlo simulation is exactly that of a K-M model in the instable 
(m^ < 3) region, with a quartic term (like the above described mean field solution) 
plus the gauge self-interaction term. Taking into account the obvious uncertainties 
of our mean field estimate the numerical agreement on the coefficient of the quartic 
contribution between our result and that of Jl^ is impressive. 



6. Conclusions 

In this letter we have shown that, besides the usual identification with a LGT 
coupled with adjoint matter in the strong coupling regime, the Kazakov Migdal 
model can also be related to the behaviour of a pure lattice gauge theory (with 
Wilson action) in the limit of very high temperature. In particular we have seen 
that the fields describe the small fiuctuations around the frozen position of the 
Polyakov loops. 

We think that this new point of view is important if we want to understand 
some of the peculiar features of the K-M model (superconfinement, lack of ordinary 
perturbative vacuum, presence of an unstable phase for < d and of a finite gap in 
the free energy), but besides this we think that there are two other reasons of inter- 
est. First the master field solution of the K-M model could give us a powerful tool 
to do both analytical and numerical calculations in pure LGT at high temperature. 
Second, more ambitious point: the K-M model is related to non critical strings. 



although this relation is still not quite clear for d > 1 (see however |]T9| and [2C 
for some recent progress on this subject). Our interpretation might lead (at least 
in the high temperature, large N limit) to a first description in terms of strings of 
LGT, a hope which justifies by itself further efforts toward a better understanding 
of the K-M model. 



Note added 

While completing this letter we received a new interesting paper by Dobroliubov, 
Kogan, Semenoff and Weiss in which the interpretation of the "filled Wilson 
loop" as the result of the perturbative expansion of gauge self-interaction term is 
discussed. The two regimes discussed in pl| may be identified within our approach 
as the high temperature deconfined regime and the low temperature confined phase 
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respectively, and the coupling A of |]21]] with our ratio |^ [|. 
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